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We studied the correlated quasi-one-dimensional systems by one-loop renormalization 
group techniques in weak coupling. In contrast to conventional g-ology approach, we for- 
mulate the theory in terms of bilinear currents and obtain all possible interaction vertices. 
Furthermore, the one-loop renormalization group equations are derived by operator prod- 
uct expansions of these currents at short length scale. It is rather remarkable that these 
coupled non-linear equations, after appropriate rescaling, can be casted into potential flows. 
The existence of what we nicknamed "RG potential" provides a natural explanation of the 
emergent symmetry enhancement in ladder systems. Further implications arisen from the 
RG potential are also discussed at the end. 

§1. Introduction 

Quasi-one-dimensional (QID) systems have attracted extensive attentions from 
both experimental and theoretical aspects.^''^ Due to the low dimensionality, strong 
quantum fluctuations often give rise to surprising behaviors, which are rather dif- 
ferent from our intuitions built in higher dimensions. It is then exciting to explore 
various exotic phenomena, such as spin-charge separation,'^''^ unconventional elec- 
tron pairingjl^'l^ spin liquids,^''^ broken time-reversal symmetr}^^'^ and so on, in 
the correlated QID systems .El El CIJ CI EHHMJMJE) 

Note that, in addi- 
tion to the ladder compounds, carbon nanotubes, nanoribbons and quantum wires, 
after integrating out fluctuations at higher energy, are also described by the QID 
theory. Therefore, not only posting a challenging task for academic curiosity, the 
understanding of the QID systems now becomes crucial as the technology advances 
to the extremely small nanometer scale. 

Since mean-field theory is not appropriate in the presence of strong quantum 
fluctuations, the ground state properties are studied by combinations of powerful 
techniques such as exact solution (Bethe Ansatz), renormalization-group (RG) anal- 
ysis, bosonization and numerical simulations. In this article, we would mainly focus 
on the RG analysis for QID correlated systems and reveal several interesting hidden 
structures behind the complicated RG flows. However, we would like to emphasize 
that RG analysis only gives rise to various instabilities of the homogenous Fermi 
liquid but does not provide a complete description of the interaction-driven phases. 
A non-perturbative complementary approach like bosonization is necessary to un- 
derstand the physical properties in low-energy limit. Since the combination of RG 
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analysis and bosonization have become rather standard now, interested readers are 
directed to several useful reviews in the literature. 

In the following, we would give a brief introduction to the current understanding 
of the ground states in QID systems. In some of the ladder compounds, the ground 
state is an (charge) insulator. Since the charge sector is frozen out, the ground-state 
phase diagram is greatly simplified and can be determined solely from the remaining 
spin degrees of freedom.!^ Early theoretical studies for Heisenberg ladders revealed 
an interesting odd-even effect .'^''^''^'^ If the number of legs is even, the 
system is expected to be a spin liquid with a finite gap to all spin excitations. On 
the other hand, if is odd, the ground state has quasi-long-range antiferromagnetic 
order and the low lying excitations are gapless and belong to the universality class 
of the single Heisenberg chain. 

The subtle interplay between charge and spin fluctuations bring out rich physics 
and exotic surprises. For instance, it has been speculated for a long time that a 
doped spin liquid may give rise to unconventional superconductivity. Let's takle the 
simplest and well-studied two-leg ladde^'^'^'-^^^^^^^ as an example. 
Without doping, there is one electron per site on average (usually referred as "half- 
filled"). Because of the Coulomb repulsion, charge excitations acquire a gap which 
makes the two-leg ladder a Mott's insulator. Besides, due to the tendency for singlet 
bond formation across the rungs of the ladder, spin- liquid behavior is expected. Both 
numerical and analytical approaches support that the ground state at half filling is 
in the Mott-insulating spin-liquid phase.'^'t^''^''^ Upon doping, the charge gap 
vanishes, these preformed pairs give rise to the quasi-long-range superconductivity 
order. It is rather interesting that the pairing symmetry of the Cooper pairs is the 
unconventional d-wave, in contract to the s-wave symmetry for phonon mediated 
superconductivity. This provides a concrete example for unconventional pairings 
and thus superconductivity in electron gas in the presence of repulsive interactions. 

Although our understanding of the strictly one-dimensional systems is greatly 
benefitted from exact solutions, it is often found that the QID systems are not 
soluble analytically. Furthermore, since the hopping along the transverse direction 
is relevant, the physics phenomena for QID systems can be dramatically different 
from the strictly ID systems. Therefore, the most reliable approach to clarify the 
competitions among various ground states is the renomalization group (RG) anal- 
ysis. Except for some exactly soluble models, the phase diagram of QID systems 
is obtained by first spotting the instabilities in RG flows and then analyzed with 
non-perturbative bosonization techniques. For A^-leg ladders, one typically needs 
to solve order A^^ coupled non-linear first-order differential equations to obtain the 
flows. Since the number of allowed interactions are large, the derivation of the flow 
equations for all couplings under RG transformations becomes formidable. On the 
other hand, both numerical and analytical approaches seem to indicate rather simple 
phase diagrams for generic QID systems. 

As one of the main points in this article, we will address the issue why the phase 
diagrams for quasi-one-dimensional systems are rather simple, while the renormal- 
ization group equations behind the scene are non-linear and messy looking. The 
puzzle is answered in two steps - we first demonstrate that the complicated cou- 
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pled flow equations are simply described by a potential V{hi), in an appropriate 
basis for the interaction couplings hi. The renormalization-group potential is explic- 
itly constructed by introducing the Majorana fermion representation. The existence 
of the potential prevents chaotic behaviors and other exotic possibilities such as 
limit cycles. Once the potential is obtained, the ultimate fate of the flows are de- 
scribed by a special set of fixed-ray solutions and the phase diagram is determined 
by Abelian bosonization. Extension to strong coupling regime and comparison with 
the Zamolodchikov c-theorem will also be discussed. 

The simplicity of the phase diagram can be partially understood by the widely 
used scaling Ansatz of the couplings in weak coupling!^ '^'^^ 

[Id - L) 

where Gi are order one constants satisfying the non-linear algebraic constraints (dis- 
cussed later) and Id is the divergent length scale where the flows enter the strong 
coupling regime. The Ansatz was motivated by the numerical observation that the ra- 
tios of renormalized couplings reach constant, as long as the bare interaction strength 
is weak enough. However, it is still puzzling why the phase diagrams, generated by 
many coupled non- linear flow equations, do not reflect the same level of complexity. 
In fact, for a complicated system with many couplings, the coupled non-linear differ- 
ential equations are likely to produce chaotic flows generically. Even if the flow is not 
chaotic, it might as well rest on limit cycles. This peculiar possibility for quantum 
systems was addressed by Wilson and collaborators in a recent paper.^ 

So the question remains: Why are the phase diagram and the RG flows so simple 
(without chaotic flow or limit cycle) in QID systems? We found the question can 
be answered by combining the weak-coupling RG analysis together with the non- 
perturbative Abelian bosonization technique. Note that the phases of the correlated 
ground state often lie in the strong coupling, so weak-coupling RG alone can not 
pin down the phase diagram. It is the powerful combination of the perturbative RG 
analysis and the non-perturbative bosonization technique which deliver the desired 
answer here. Note that, in the low-energy limit, the QID systems involve Nj fla- 
vors of interacting fermions with different Fermi velocities, where Nf is the number 
of conducting bands. Although constrained by various symmetries, the number of 
allowed interactions is tremendously large as A'^^ increases. In general, the RG equa- 
tions to the lowest order are already very complicated - not to mention solving them 
analytically. 

However, quite to our surprises, we found that, at one-loop level, the RG flows 
can be derived from a potential, i.e. the coupled non-linear flow equations can be 
cast into this elegant form by an appropriate choice of coupling basis, 

dhi _ dV{hj) 
dl dhi 

where V{hj) is the RG potential. We emphasize that this is only possible after a 
unique transformation of the couplings, hi(l) = LijQj^l) (up to a trivial overall factor 
for all couplings), where Lij is some constant matrix. The existence of the potential. 
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which requires the coefficients in the RG equations to satisfy special constraints, also 
provides a self-consistent check on the RG equations derived by other approaches. 

The flows of Eq. Hl-2|) in the multi-dimensional coupling space can be viewed as 
the trajectories of a strongly overdamped particle moving in a conservative potential 
V{hi). Note that the change of the potential V{hi) along the trajectory is always 
decreasing, 



that the function V{hi) never increases along the trajectory and is only stationary 
at the fixed points where dhi/dl = 0. Thus, the RG flows have a simple geometric 
interpretation as the trajectory of an overdamped particle searching for potential 
minimum in the multi-dimensional coupling space. 

This simple geometric picture rules out the possibilities of chaos and the exotic 
limit cycles in QID systems. The ultimate fate of the flows would either rest on the 
fixed points or follow along the "valleys/ridges" of the potential profile^^' to strong 
coupling. Since there is only one trivial fixed point at one-loop order, most of the 
time, the flows run away from the non-interacting fixed point. Starting from weak 
enough bare couplings, the ultimate fate of the flows is dictated by the asymptotes 
of the "valleys/ridges" of the potential profile. It provides the natural explanation 
why the ratios of the renormalized couplings reach constant in numerics. That is to 
say, the existence of RG potential implies that the ultimate fate of RG flows must 
take the scaling form described in Eq. Hl-ip . Detail properties of these asymptotes, 
referred as fixed rays, will be discussed in later section. 

Since the ultimate fate of RG flows is described by the simple Ansatz in Eq. , 
the specific ratios of couplings simplify the effective Hamiltonian a lot. Making 
use of the Abelian bosonization, one can determine which sector acquires a gap, 
triggered from the weak-coupling instability. The phase of the ground state is then 
determined by watching which fixed ray (asymptote) the flows go closer to. Because 
there are only limited solutions of the fixed rays, the phase diagram is rather simple. 
Therefore, by combining the powerful techniques of weak-coupling RG and Abelian 
bosonization, we pin down the reason behind the simple-looking phase diagram out 
of the messy non-linear fiow equations. 

In fact, the combination of weak-coupling RG and Abelian bosonization goes 
beyond the usual mean-field analysis and is crucially important when there are more 
than one competing orders.!^ For instance, Lauchli, Honerkamp and Rice recently 
studied the so-called "d-Mott" phase in one and two dimensions, where antiferromag- 
netic, stagger-flux and d-wave pairing fluctuations compete with each other simul- 
taneously. The conclusion drawn from the numerical density-matrix RG in strong 
coupling agrees rather well with predictions made from one-loop analysis in weak 
coupling. This lends support to the powerful combination of weak-coupling RG and 
Abelian bosonization approach for strongly correlated systems. Since the method 
of bosonizing the flxed rays is already developed in previous papers we 
would concentrate on the novel existence of RG potential in this paper. In particular. 
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we would construct the RG potential explicitly. 

The existence of RG potential and its elegant geometric interpretation are also 
related to the mysterious emergence of symmetry enhancement in ladder systems in 
low-energy limit. Let us return to the simple example of the half- filled two-leg ladder 
again. In general, one expects that charge and spin gaps could be rather different. 
Surprisingly, a complete degeneracy of charge, spin and single particle gaps emerges 
in the asymptotic weak coupling limit. Based on the one- loop renormalization group 
(RG) analysis, the effective low-energy theory for the two-leg ladder is described by 
the exactly soluble Gross- Neveu model with a global SO (8) symmetry. Other than 
the exact degeneracy of excitation gaps, many interesting results can be drawn from 
the Bethe Ansatz solution. It turns out the emergence of the enhanced symmetry 
is directly tied up with the "valleys/ridges" of the RG potential profile in two-leg 
ladder. One can show analytically that these "valleys/ridges" in potential profile 
correspond to asymptotic flows with fixed ratios between all couplings that give rise 
to the enhanced symmetry. 

There have been criticJ^'^ about the unexpected emergence of the enhanced 
SO (8) symmetry, questioning the stability of the symmetry under generic perturba- 
tions and the limitation of the RG equations derived in weak coupling. It was argued 
that the emergent symmetry derived by perturbative calculations is fragile because 
the system eventually flows toward the fixed point in strong coupling. Both con- 
cerns can be addressed by performing a complete analytical study for the stability 
in the vicinity of the symmetric phases. It should be emphasized that this emergent 
symmetry enhancement in weak coupling does not implies that the fixed point in 
strong coupling is S0(8) symmetric at all. On the contrary, the symmetry reflects 
the local topology of RG flows near the trivial Fermi liquid fixed point which can 
be safely described by perturbative calculations. Indeed it is rather obvious that, 
when the interaction U is much larger than the bandwidth t, the charge gap is much 
larger than the spin gap. Thus, the symmetry in the strong coupling is completely 
destroyed. As to the stability of the symmetry under generic perturbations, the 
complete stability check shows that the SO (8) symmetry is robust in weak coupling. 
Even though some couplings are relevant according to the conventional classification, 
they do not destroy the symmetry but only give rise to anomalous corrections which 
scale as [U/ty/^. A new classification of "relevant" and "irrelevant" perturbations 
is necessary here. 

The rest of the article is organized as following. In Sec. 2, we present the model 
for correlated quasi- ID systems and derive the effective theory in weak coupling. 
Making full use of SU(2)xU(l) symmetries, all interaction vertices can be elegantly 
expressed in terms of various bilinear currents. In Sec. 3, the one-loop RG equa- 
tions for the correlated quasi-lD systems are derived from the technique of operator 
product expansions of bilinear currents. In Sec. 4, we used the half-filled two-leg 
ladder as example to reveal a hidden structure behind the messy RG equations. By 
appropriate rescaling of all coupling constants, we show that the RG fiows can be 
derived from a single potential. We also demonstrate that the RG potential also 
exists for the doped ladders. In Sec. 5, we provide a formal proof for the existence of 
the RG potential by transforming into the so-called Majorona basis. The potential 
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gives an elegant geometric interpretation of destinations of the RG flows. Finally, 
in Sec. 6, we address the interesting issue about the emergence of enhanced sym- 
metry in ladder systems and its connection to the existence of RG potential. We 
also present a simple scaling argument near these symmetric rays and argue that 
anomalous scaling is expected in the vicinity. 



§2. Quasi-ID Ladders 



We start with the simplest QID ladders. In tight-binding limit, the Hamiltonian 
of the N-\eg ladder, H = Hq+Hj, contains the nearest-neighbor hopping and generic 
short-range interactions. 



Ho=Yl\- Mai^ + ^)dzaix) + dl^{x)dia{x + 1)] 
x,i,a 

-t±[dl+^,^{x)di^{x) + (2-1) 
Hi= Yl niix)Vi,j{x-x')nj{y), (2-2) 



X,X' ,l,J 



where di{d-) is the annihilation (creation) operator for the fermion on chain i (i = 
1...N) and a =],[ is the spin index. The parameters t and t± denote the hopping 
amplitudes along the same chain and across the rungs between neighboring chains. 
The density operator on the ladder is ni{x) = ^a'^la(-^)'^io!{x)- The short-range 
interaction Vij{x — x') is taken to be the density-density type here but it can be easily 
generalized. Note that the method we try to build here is for generic interaction 
profile, not limited to the on-site repulsion. 

In weak coupling, it is natural to concentrate on the hopping Hamiltonian Hq 
first, before including the correlation effects from Hj perturbatively. Depending on 
the boundary condition in the transverse direction, the nearest-neighbor hopping can 
be diagonalized easily. Since the open boundary condition (OBC) is more natural 
and relevant to the realistic ladder materials, we will concentrate on this case first. 
However, we also include the studies on periodic boundary condition (PBC) later 
for completeness. For OBC, the transverse hopping across the rungs can be brought 
into diagonal form by the N x N matrix Sjm, 



djaix) = Sjmil^rna{x), with Sjm = J sin i j^—j j . (2-3) 

Furthermore, after transforming into the momentum space along the chain direction, 
the hopping Hamiltonian takes the simple form, 

^0 = E r 1^ ^^ik)^Lik)i^^a{k). (2-4) 

The single-particle spectrum is ei{k) = —2t cos k — 2t± cos{kyi), where the quantized 
momentum in transverse direction is kyi = m/^N + 1). 
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For the purpose of deriving effective theory in low-energy hmit, it is sufficient to 
linearize the energy spectrum near the Fermi points that intersect with the chemical 
potential ei(fc_pj) = fi. This implies the lattice operator ipi^ix) can be decomposed 
into a pair of chiral fields at each Fermi points, 



(2-5) 



Since the fast-oscillatory phase factors near the Fermi points are pulled out already, 
the chiral fields iljpi„{x) with P = R,L are smooth varying operators. Finally, the 
hopping Hamiltonian can be expressed in terms of these chiral fields with linear 
dispersion, 



Ho 



dx Vj 



(2-6) 



where Vi = 2tsm{kFi) is the Fermi velocity of the i-th band. 

Before we dive into the details of how to derive all possible interactions, we 
would like to comment on the validity of the weak-coupling theory first. To address 
this issue, let's write down the partition function for the correlated QID ladder first, 

Z =Tr e-l^" = J L»[V^V]e"'^. 
The Euclidean action can be derived from the path integral, 



/ dT 

Jo 



x,i,a 



(2-7) 



The validity of chiral-ficld decomposition relies on how the interaction strength is 
renormalized after integrating out the modes that are far from the Fermi points. To 
be more precise, the momentum regimes \ p — kpi \> A and \ p + kpi \> A that are 
outside the vicinity of the Fermi points, are integrated out in the path integral. The 
standard perturbation theory gives the renormalized interaction strength 



Urc^U 



1 -I- const X — m — — 

t \ A 



(2-8) 



Therefore, as long as {U/t) <€i l/ln(^) ~ 0(1), the one-step RG from the whole 
Brillouin zone down to the cutoff A does not give significant renormalization and the 
perturbative approach is reliable. 

We are now ready to search for all possible interaction vertices and write them 
down in terms of the chiral fields. Note that, at half filling, there are special con- 
straints between Fermi momenta, kpi + kp^ = vr, where i = N + 1 — i. These con- 
straints give rise to additional vertices known as umklapp interactions. In addition 
to the U(l) X SU(2) symmetry corresponding to charge and spin conservation, these 
terms must be preserved by charge conjugation, time reversal, parity, and (lattice) 
translations. The most general particle-conserving 4-point vertex takes the form. 



(2-9) 
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Here the spin indices are suppressed for notational clarity. We will include them 
back explicitly later. 

Various discrete symmetries, such as charge conjugation, time reversal, parity 
and so on, give constraints on Pa,ia- However, the strongest limitation arises from 
conservation of lattice momentum, imposing constraints on the allowed band indices 

{Pa,ia) 

- Pikpi, - P2kFn + P^^Fvi + Pikpu = 2^n^. (2-10) 

Although we deal with OBC here, it is important to remind the readers that the 
transverse momentum is also conserved when PBC is considered. Thus, there are 
additional constraints associated with the translational invariance in the transverse 
direction, 

- Pih - P2i2 + Pzh + PiH = Nuy (PBC only!!). (2-11) 

The remaining work is to find out all possible solutions for {Pa,ia) in Eqs. (|2-lUp - 
H2-ll() . When searching for the solutions, one must keep in mind that a solution 
with the same chirality for all Pa does not cause any instability in the ground state. 
Instead, it only gives rise to renormalization of Fermi velocities and, therefore, can 
be safely ignored for the derivation of one-loop RG equations. Since they are irrel- 
evant to the phase diagram of ground states, we will ignore them in the following 
discussions. 

2.1. Open Boundary Conditions 

For OBC, only momenta in the kx direction is conserved. Therfore, only Eq. (|2-in|) 
needs to be studied. To pin down all possible solution is literally impossible. How- 
ever, one notice that we can classify the solutions into two categories: the major 
vertices that always exist for generic parameters t,t±,n,... and the minor vertices 
that only show up when the parameters are specifically fine-tuned. Since the phase 
space of the minor vertices are negligible in comparison with the major ones, it is 
reasonable to drop them. However, one should keep in mind that the approxima- 
tion may break down if the system under consideration has peculiar shape of Fermi 
surface that greatly enhances the minor vertices. The reason why the major ver- 
tices appear in the whole parameter space implies the momentum conservation in 
Eq. I|2-1U() must be satisfied trivially. Thus, we look for linear combinations of some 
trivial identities for half-filled QID ladders. For OBC, these identities are rather 
simple, 

kpi - kpi = 0, (2-12) 
kpi + k J,. =Tr, (2-13) 

where i = {N -|- 1) — i is the band index which has the same Fermi velocity as the 
i-th band. The first identity is indeed trivial and the second one is protected by the 
particle-hole symmetry at half-filling. 

Vertices can be classified by the integer in Eq. (|2-1U() . Since kpi < tt, Ux can 
take values, 0, ±1. For = 0, we can construct the linear combination of Eq. (|2-12|) . 
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kpi — kpi + kpj — kpj = 0. (2'14) 

Compared with Eq. H2-10p . it is straightforward to spot two sets of solutions. The 
first set comprises forward scattering that satisfies 

(Pl,P2) = (P3,A) 

(11,12) = {is, h), (forward), (2-15) 

The second set is the Cooper (or backward) scattering, 

Pi = P2, -P3 = Pa', 
n = ^2, h = U, (Cooper). (2T6) 

In Eqs. H2-15() - H2-16|) . and in the remainder of the paper, P = —P and {xi,X2) = 
(x3,X4) indicates pairwise equahty, i.e. either xi = X3, X2 = X4, or xi = X/i, X2 = x^. 
The two possible solutions for forward scattering actually describe the same vertices, 
up to a sign from the fermion ordering. Referring to the 2D Brillouin zone, one sees 
that forward scattering preserves the particle number in each band separately, i.e. 
electrons are annihilated/created within the same band. For Cooper scattering, 
pairs of electrons are annihilated in one band (for instance, the ii-th band) and then 
scattered into anotehr band (for instance the is-th band). 

There is another set of solution for = 0. Making use of the linear combination 
of Eq. (imijl 

kpi + kpi-kFj -kp~. = Q, (2-17) 

we found another set of vertices which involve /our different bands. In contrast to the 
forward and Cooper scatterings that only involve two different bands, this strange 
set of solutions requires 

(Pi, P2) = (P3, A), Pi = P2, Pa = P4; 
{ii,i2) = (iA,n)- (2-18) 

Note that the existence of this strange set of vertices crucially relies on the nested 
Fermi surface at half filling. However, one should not confuse it with the umklapp 
interactions (where rix 7^ 0) because the momenta in these strange vertices are 
exactly conserved. 

In addition to the three set of vertices for rix = We now turn our attentions 
to the so-called umklapp interactions with non- vanishing rix- For rix = ±1, the only 
useful combination from Eq. ()2-13() is, 

kpt + kp:. + kpj + kp~. = 27r. (2-19) 

Compared with Eq. H2-lUp . it gives two sets of solutions. The first set is the usual 
umklapp interaction that requires 



Pi = P2 = Pi = Pa; 
{11,12) = {h,H)- 



(2-20) 
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This set of vertices describes the scattering process for a pair of electrons from 
bands i and j with one chirahty to their particle- hole symmetric partners, bands i 
and j, with opposite chirality. The change of momentum associated with this kind 
of vertices is not zero, AP = ziz(kFi + kpj + kp~- + kp-j) = ±2?! and thus belong to 
the category of umklapp interactions. 
The other set of solutions requires 

Pi = P2 = P3 = Pa; 
h = ^2 h = U- (2-21) 

This set of vertices describes the scattering process for a pair of electrons from 
particle-hole symmetric bands i and i with one chirality to another pair of bands j 
and j with opposite chirality. Again, this kind of vertices falls into the category of 
umklapp interaction because the change of momentum associated with this kind of 
vertices, AP = ^{kpi + kp-- + kpj + kp-j) = ±2it, is not zero. 

Collect all the results we got so far, the allowed interactions consist of three 
kinds of momentum-conserving vertices and two kinds of umklapp interactions. It is 
possible to write down the Hamiltonian density Hint = + Ti-l^^ in terms of the 
chiral fields, 

n^t = E {c[P,i,MUi^pA^^Tj 

P,i,j ^ 

+ F[p,i,i]V'k^«4iV'Pi 



P,i,j ^ 



+ W[P,^J]i;UiJl,^PJ,.^P-^. (2-23) 



Note that Wj-^) contains vertices which conserve momentum, while H^-^^ denotes the 
umklapp interactions in the kx direction. 

Now we need to include the spin indices. If we follow the conventional g-ology 
approach, it is rather messy to write down all vertices explicitly. To make the full 
use of the rotational symmetry in spin sector, it is very helpful to express the four- 
fermion interactions in terms of SU(2) bilinear currents. We introduce the SU(2) 
scalar and vector currents, 

Jpij = ^'^Uc.'4'pjc, •^p«j = ^ V'k^o-a/jV'pj^), (2-24) 

^Pij = ^'^Piceal3lf^Pji3, Ipij = ^ i^Pic{eO-)al3lpPj/3, (2-25) 

where cr is Pauli matrices and e is the Levi-Civita antisymmetric tensor with the con- 
vention ei2 = —621 = 1- To regularize the composite operators in Eqs. ()2-24|l - ()2-25|) . 
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the defined currents are already normal ordered (although we do not indicate the nor- 
mal ordering explicitly). As is clear, Jjj, lij are scalars under SU(2) transformation, 
while Jij, lij are vector currents because of the Pauli matrices cr. Furthermore, due 
to Fermi statistics, it is straightforwad to show the scalar current lij is symmetrical, 
while the vector one lij is anti-symmetrical, 

lij — Iji, lij — ^ji- (2'26) 

Therefore the diagonal part of the vector current In is identically zero because of 
the Fermi statistics. 

The four-fermion interactions can be expressed in terms of products of the bilin- 
ear currents. Since the Hamiltonian is an SU(2) scalar, only products of scalar-scalar 
currents or vector-vector ones are allowed. Putting spin indices back into the theory, 
the momentum-conserving vertices in Eq. ()2-22|) can be classified into spin/charge 
and scalar /vector sectors, 

^(1) ^ jR jL _p<T jR . jL 

' ^int ij ij U «J V 
_l_ fp jR jL _ la jR _ jL 

J ij a jj Jij'^ a ' jj 

+ ^ij-^ij'^]i ~ ^ij-^ij ■ "^fi' (2-27) 

where we have used Einstein's convention that summations over the repeated indices 
are implied. The first two sets of couplings fij and Cij denote the forward and Cooper 
scattering amplitudes between bands i and j at generic fillings. However, at half 
filling n = 1, due to the particle- hole symmetry near the Fermi surface, additional 
vertices Sij are also allowed. Note that the vertices Sij involve four different bands 
and are rather different from the forward and Cooper scattering, fij and Cij, at 
generic fillings which only involve two bands. 

Although the vertices that conserve momentum can be elegantly written down 
explicitly in terms of SU(2) currents, not all of the couplings are independent. Some 
of the vertices are doubly counted and should be removed. Besides, various discrete 
symmetries would further reduce the number of independent couplings. This is 
actually a good news since we eventually need to solve these couplings numerically. 
To avoid double counting, one notices that fnjCn describe the same vertex. Thus, 
we choose the diagonal elements of the forward scattering to be zero, i.e. fa = 0, 
to avoid double countings. For the same reason, we choose to set sn = = s---. 
Furthermore, under charge conjugation, the indices in the currents switch Jij Jji, 
which implies Cij = Cji and Sij = Sji. Following similar arguments, parity symmetry 
(in x) implies fij = fji. Combining the above symmetrical properties with particle- 
hole symmetry, we found that Aij = A--^ where A = {c, /, s}. While it is not obvious 
at this point, the choice of signs for the scalar and vector couplings in Eq. (|2-27() is 
such that they are all positive for the repulsive on-site interaction. 

Now we turn to the umklapp interactions which do not conserve momentum 
exactly (but up to reciprocal lattice vectors). Comparing the difference between 
interactions in Eqs. ()2-22|) and ()2-23|) . we need the other set of SU(2) currents Ipij 
and Ipij to construct all possible umklapp interactions. Again, the allowed four- 
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fermion interactions can be expressed in terms of the products of SU(2) currents, 



Since ua^wu describe the same vertex, we choose wa = = w--- to avoid double 
counting. Besides, the product of currents are symmetrical (even though single lij 
is anti-symmetrical), we can choose the couplings to be symmetrical Uij = Uji and 
Wij = w--'a- Under parity transformation and charge conjugation, it is easy to show 
that Uij = U--J and Wij = Wji. From Fermi statistics, it is easy to show that u^- = 
and = 0, where m = (iV + l)/2 for OBC's and m = N/i for PBC's (although 
we haven't discussed the details yet). So far, we only concentrate on the OBC's. 
In principle, the change of boundary conditions can lead to different sets of allowed 
interactions in ladder systems. However, as we will explain in the following, it is 
rather surprising that the above interaction H = H^l^^ + wj^l derived for OBC's is 
also correct for PBC's. 

2.2. Periodic Boundary Conditions 

For PBCs, the interactions strongly depend on whether the number of chains is 
odd or even. For odd-chain system, half filling is not even an umklapp line. As a 
result, the allowed interactions are greatly simplified and only Cooper and forward 
scattering in Eq. (|2-27() are allowed. The detail ground-state phase diagram and RG 
analysis can be found in the literature.^ When the number of chains is even, the 
situation is much more messy. The complication arises from more identities at half 
filling. 



where i = sign(i){N/2) — i. The it sign in the identities and the dependence of sign(i) 
in the definition of i makes the analysis a bit messy. Note that these additional 
identities reflect the peculiar properties of the nested Fermi surface at half filling. 

As in previous subsection, we classify the vertices by {nx,ny), where rix , riy can 
take values, 0,ibl,ib2. However, Ux = ±2 is ruled out because it corresponds to 
the interactions between completely filled bands which, as mentioned before, do not 
participate the low-energy physics here. For Uy = ±2, all the bands are lying on 
the boundary of Brillouin zone, i.e. ia — which actually indicates the same 

band and thus included in the Uy = case. Therefore, only four kinds of vertices 
are allowed: (0,0), (0,±1), (±1,0), and (±1,±1). 

For {ux^Uy) = (0,0), the momenta is conserved in both directions. We can 
construct the linear combination from Eq. ()2-29() 




(2-28) 



kpi - kp±i = 0, 
kpi + kp^i = vr. 



(2-29) 
(2-30) 



kpi - kF±i + kpj - kF±j = 0, 



(2-31) 



RG Potential for Quaisi-ID Correlated Systems 



13 



which gives rise to two sets of solutions. The first set of the solutions is, 

{ii,i2) = {±i3,±H), (2-32) 

and the second set is, 

Pi = P2, P3 = Pi', 

ii = ±12, h = ±H- (2-33) 

However, since the transverse momenta is also conserved, the solutions need to satisfy 
Eq. (|2-11|1 as well, which fixes the ambiguity of the signs. The resulted vertices are 
again the Cooper and forward scattering. 

Following the same spirit, other interactions can be constructed by different 
linear combinations of the identities arisen from the nested Fermi surface. The 
linear combination from Eq. ()2-30() gives 

kpi + kpj^i - kpj - kp^j = 0. (2-34) 

The conservation of transverse momenta puts a more severe constraint in this case. 
For forward and Cooper vertices, it rules out the annoying sign ambiguity and leads 
to the same vertices as in OBC's. Here, not only pin down the signs, it also puts 
constraints on the band indices 

Pl=P2, P3 = P4, 
(Pl,P2) = (P3,P4), 

(n,Z2) = (14,^3), ii-Z2>0, (2-35) 

where 11,12 must have the same sign! This strange constraint is quite tough to handle 
when one wants to perform RG analysis. But, as will become evident in a moment, 
The other missing vertices with opposite signs ii ■ 12 < appears in another sector 
and the strange constraint is indeed an artifact of how we classify the vertices. 

It is quite similar to find the (0, ±1) vertices. In stead of repeating the whole 
calculations again, we briefly sketch the steps to obtain the solutions. First, the 
conservation of kx momenta is the same as for (0,0) vertices. The difference arises 
because the right-hand side of Eq. ()2-ll|) is not zero. After a liitle algebra, three sets 
of solutions are found. The first set of solutions is 

(Pl,P2) = (^'3,^4), 

- - TV 

{h,i2) = \h - ^21 = —; (2-36) 

and the second one is. 

Pi = P2, P3 = Pi] 

N 

k=i2, i3 = k, Ih - i3\ = TT ■ (2-37) 
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The last one is 



Pi = P2, 

iPl,P2) 

{h,i2) 




h- 12 < 



(2-38) 



where ii,i2 must have the different signs. One notices the last set of solution has 
appeared in the (0,0) vertex except the constraint on the band indices is different. 

It is straightforward to show that the first two sets of solutions are actually 
included in the last one already. With the simplification, we can write down the 
Hamiltonian density for (0, 0) and (0, ±1) vertices together in terms of products of 
SU(2) currents. Now it is clear that these vertices are described by 7^1^) in Eq. (|2-27l) 
- the same as in OBC's!! Similarly, when combining the (±1,0) and (±1, ±1) vertices, 
we end up with the Hamiltonian Til^^ in Eq. (|2-28() . It is worth commenting that 
the extra subtleties coming from the conservation of transverse momentum make the 
search for all allowed vertices much more difficult and messy. But, at the end, the 
effective Hamiltonian takes exactly the same form as in OBC's, except the definition 
of dual points on the Fermi surface are different. At the time of writing, we do not 
fully understand the reason behind this coincident. 



In stead of following the conventional diagrammatic technique, we derive the one- 
loop RG equations by operator product expansions (OPE) of the bilinear currents. 
Since the leading divergence is logarithmic, the one-loop RG equations are universal 
despite of different schemes. In fact, the universality is crucial since we later show 
that the exponents of anomalous scaling are directly determined by the coefficients 
in the RG equations. Since the exponents of various power laws are robust, we 
expect the coefficients of the one-loop RG equations to be universal, which agrees 
with the logarithmic divergences of marginal interaction (four-fermion vertices) in 
QID systems. 

In the following, we outline the steps to obtain the RG equations and also derive 
one of the renormalized coupling as a simple example to familiarize readers with the 
technique. To proceed, it is convenient to introduce the partition function in path 
integral formalism. 



where Zq is the partition function without interactions and (e~'^*"')o with the sub- 
script denotes the average over the unperturbed action Sq only. We can reexpo- 
nentiated it by the cumulant expansion, 



§3. One-Loop RG 



Z = 
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The only nontrivial term is {Sf^^)o which will renormalize the other interactions. To 
one-loop order, 

/•oo ^ 

{-Sfnt)o = / Yi dXi {'Hint{Tl,Xi)'Hint{T2,X2))Q , (3-3) 

i=l,2 

the notation = Ja<\xi-'X2\<B' ^^'^^ integral has two parts: long-wavelength 
modes, 6a <| xi — X2 |< oo, and short- wavelength modes, a <\ xi — X2 \< ba, where 
6 > 1 is the rescaling parameter. When we integrate out those short-wavelength 
modes, since all fields are at nearby space points, wc can use the OPE to replace the 
product of them by a series of local operators, that is, 

/ba p 
(iTi dXi {nintiTl,Xi) nint{T2,X2) )o - dv dx {6nint)o ■ (3-4) 
- i=l,2 

The effective interaction SHint has the same form as the original Hamiltonian and 

renormalize corresponding couplings when the short-range fluctuations are progres- 
sively integrated out. Studying how the couplings get renormalized, the one-loop 
RG equations can be derived. 

Here we demonstrate how to calculate the effective Hamiltonian SHint with the 
use of OPE's explicitly. First of all, when chiral fermion fields ipR/Lia (with the same 
chirality) are brought closer together in space and time, the leading singularities are 
simple poles, 

V'ffia(x,r)Vle,-^(0,0) ~ ^ + 0(1), (3-5) 

Zi 

i^Uaix, r)4^.^(0, 0) ~ ^ + 0(1), (3-6) 

where zi = 2n (viT — ix). The above equations are valid as long as (x, r) close enough 
to (0,0). The full set of OPE can be obtained, making use of the above identities. 
For instance, when two right-moving currents Jjnj and Jjiki are brought closer in 
space and time, their product can be approximated by another operator, 

Jijix,T)JkiiO,0) = v4(x,r)V',-a(x,r) ^1^(0, 0)V'i;3(0, 0) 

- -^-^4,mja{z,) + ^vL(^o^/Mo) = - ^j., + ^j..(3-7) 

Zi Zj Zi Zj 

Here, for notational convenience, the subscript of the chirality is suppressed. The 
others operator products can be calculated in the similar way. The full set of OPE 
can be found in Appendix A. 

With these OPE's, we now perform explicit calculations for a typical term in 
{-Sint)Q^ say 

UiJ^kl I {j!',{z)Jt,{z)4i{w)jj:^{w)) (3-8) 
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where ^ denotes a four-dimensional integral over two complex planes z and w. 
In the RG procedure, we integrate out degrees of freedom at short-length scale, 
i.e. where the complex varibles z and w are close together. Replacing the operator 
products by equivalent local operators, the OPE's give 



Ojk 



{Zj - Wj) 



7« 



{Zi - Wi) 



Ojk 



1 



1 



1 



J, 



kj 



8 I z, 



1 



1 



{Zi-Wi){z* -w*) 



1 



The integration over the fast modes is elementary, 

1 In 6 



f dx I 

J a<\x\<ba J ~i 



dr 



dl 



ZiZ* 7r{vi + Vj) Tr{vi + vj) 



(3-9) 



After integrating over the short-range fluctuations and collecting terms together, the 
renormalized action is, 



dl 



4 ''^ 27rvk 



1 (~p x2 dl 



(3-10) 



Compared with the bare Hamiltonian, it contributes to renormalization of Cooper 
scattering cf^- and forward scattering f^-., 

dl 



P 7.P 

4 ^ 2TTVk ' 



dl 



TT{Vi+Vj)' 



(3-11) 
(3-12) 



To simplify the notation, we define Aij = Aij/-K{vi + Vj), where A = {b, f, s,u,w} 
and Qjj^fc = {vi + Vk){vj + Vk)/2vk{vi + Vj). The RG equations are. 



p p 



(3-13) 
(3-14) 



where A = dA/dl. One can follow similar steps to compute the renormalization of 
all couplings to one-loop order. After lengthy algebra, the whole set of coupled RG 
equations is derived and presented in Appendix B. 

Note that the renormalization of interacting vertices is described by the coupled 
nonlinear differential equations such as Eqs. H3-13p and ()3-14() . From mathematical 
viewpoint, the flows of these nonlinear equations should exhibit chaotic behavior. 
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However, numerical studies of these messy equations give rather simple and robust 
phase diagrams. Somehow, the flows is "protected" from the generic chaotic fate. It 
turns out that there exists hidden structure behind the messy-looking RG equations. 
After appropriate rescaling, the whole set of equations can be derived from a single 
potential which we nickname as RG potential. As the old saying states - seeing is 
believing, we will present two examples of RG potential in the next section first. 
Then, in following up section, we give the general proof for the existence of RG 
potential but construct it explicitly. 

§4. Hidden Structure in RG Flows 

The first example wc plan to present is the two-leg ladder at half filling. After 
chiral field decomposition, the kinetic energy consists of four pairs of Dirac fermions, 

'Ho = tpl.io.{-ividx)^Ri„ + ^L(«^^i^x)V'«a, (4-1) 
where i = 1,2 denotes the band index and a =|, | are the spin indices. Due to the 
symmetry between bonding and antibonding bands, the Fermi velocities are equal 
vi = V2 = V. Since there are only two interacting bands, the interactions are not as 
complicated as we derived before, 

'Hint — ^jJmjJLij ^ij*^ Rij ' Lij + /[jJRuJLjj fij"^ Rii ' Ljj 

+ -^{iRijlLji + lujlRji) ^{^Rij ■ ^Lji + I Lij ■ Irh)- (4-2) 

Summation over the band indices i,j = 1,2 are implied. Note that the other two 
kinds of vertices, s^j'^ and w'^f , do not show up because of the number of interacting 
bands is only two here. 

Since the Fermi velocities in bonding and antibonding bands are the same, we 
can choose the Fermi velocity 2ttv = 1 to simplify the numerical factors. To give the 
readers a fiavor of the complexity of the coupled non-linear RG equations, we write 
them down explicitly here, 

^ = -{{^2? - \{cl.f + \{u{2f + \{ul,)\ (4-3) 

J^cr 1111 

= -(cTl)' - 2^2C?2 - ^{^2? - ^<2^^U - 2(^2)', (4-4) 

dc^2 1 3 1 ~ 3 ~ 

= ~2^11^12 ~ 2^11^12 + 2^2/12 + 2^12/f2 + '"'12^11' (4"5) 

dcl^ 1 1 1 ~ 1 ~ _ 

= ~ 2^11^12 ~ 2^12^11 ~ -I- 2/12^12 + 2^12/^2 ~ C12/12 + ^ii'^i2) (4'6) 

% = {{^2? + \{cl2f + \{<2f + l{ui2? + (4-7) 

^ = -ifuf + ^C?2Cl2 - l{c1,r + in?2<2 - ^«2)', (4-8) 
dil'' ~ 1 S 

^ = /f2< + lf.2u{2 + |c?2n?2, (4-9) 
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(4-10) 



1 



1 



11 "12 



,'-11 "12 



1 



1 ^ 



C11UI2 + 7T/f2^?2 + o/l2<2 - /r2<2 + C?2<- (4-11 



Neither do the above equation look inspiring, nor do they seem to reveal any struc- 
ture/symmetry behind. However, if we rescale the couphngs in the foUowing way. 



'"ll' '^12' ^1' ^2' /l2' '^12' ^11' ^12' fii 



"-11 „P J> rfi fP Ci2 fi2 

/pr) "'12' ''11' '"12' J 12' ^1 ini ini in 



LV2 



a/3' V3'x/3' Vs. 



(original gi) 



(rescaled hi). 



(4-12) 



In terms of the rescaled couplings /li, the coupled non-linear RG equations can be 
mapped into potential flows, 



dhi 
'df 



dV_ 



(4-13) 



where all the RG flows are captured by a single potential. 



1 



'12J 



1 ^ «/i n r , n -.o , ^ ^o^ 1 



+ 4(/f2 - [(C?2)' + K2)'] + jif?2 + C?l) [K2)' + K2)'] 

+ \ifl2 - CW [CIA2 + <2<2] - ^ifl2 + cTl) [(C^s)' + «2)'] 
+ -^<,<2C?2 + ^<1<2C?2 + ^Kl)Vf2- 



(4-14) 



This is truly remarkable that the complexity of the RG equations can be removed 
by simple rescaling. Since the flows under the influence of potential can be viewed 
as the trajectories of an overdamped particle in multi-dimensional coupling space, 
we know the ultimate fate of the RG flows is nothing but chasing the minima of the 
potential V . To one-loop order, the potential does not have any minimum (except 
the trivial one at hi = 0). Therefore, the flows are attracted to the fastest descending 
asymptotes that eventually go into the strong coupling regime. It is rather easy to 
find out the asymptotes where the couplings keep constant ratios. That's why chaos 
are never spotted in numerical studies for the RG equations of QID systems. 

One may question that the two-leg ladder at half filling may be a special case 
and the generic QID systems are not necessarily described by the RG potential. 
Therefore, it is insightful to consider the doped A'^-chain ladder at generic fillings. 
The kinetic energy is described by 2N (the factor of 2 comes from spin) pairs of 
Dirac fermions with different velocities Vi. At generic fillings, the shape of Fermi 
surface does not have any special symmetry. As a result, only forward and Cooper 
vertices are allowed, 



J 



(4-15) 



RG Potential for Quaisi-ID Correlated Systems 



19 



where summations over the band indices i,j = 1, 2, A'^ are implicitly indicated. To 
avoid repetition, we do not intend to write down the RG equations and redirect the 
interested readers to Appendix B. The RG equations for the doped QID ladder can 
be obtained by setting the other couplings s^j'^ = u^j'^ = w^'^ = in the equations. 

The presence of different Fermi velocities makes the search for appropriate rescal- 
ing factors not a trivial task, 



4 = 27r (4^2- 
cf. = 27r(4V^^)4, 



c% = 27r 



4^/2 



/5 = 27r{^V^)/5. 



(4-16) 



It is rather remarkable that, despite of the presence of difference Fermi velocities, 
the RG potential still exists! The whole set of complicated RG equations can be 
derived from the potential below 



V(9i 



i<j 
i<j 



i<j 



\/2 \/2 2 

V3 " V3 " Vs^'''^'' 



-V2cS-V2c^. + 27,,/- 



i<j<k i<j<k 



(4-17) 



i<j<k 



where the Fermi velocities enter the RG potential only through the ratio of geometric 
and algebraic averages for different pairs of Fermi velocities. 



7ij 



Vi + Vj 



(4-18) 



Since both the half-filled two-leg ladder and the doped iV-chain system at generic 
fillings are governed by the RG potential, it motivated us to look for a more general 
way to construct the potential in next section. 
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§5. RG Potential 

To prove the existence of the RG potential, it is helpful to study the general 
feature of one-loop RG equations of the QID systems first. In weak coupling, the 
most relevant interactions are the marginal four-fermion interactions, described by 
a set of dimensionless couplings gi. The RG transformation to the one- loop order is 
described by a set of coupled non-linear first-order differential equations 

^ = Mt9m^F,, (5-1) 

where the coefficients M^^ = M^-' are symmetrical by construction. These constant 
tensors M/*^ completely determine the RG flows. 

The solution for Eq. (|5T|) can be viewed as the trajectory of a strongly over- 
damped particle under the influence of the external force Fi in the multi-dimensional 
coupling space. Naively, one might rush to the conclusion that the conditions for 
the existence of a potential requires are, 

which implies that the tensor Mf^is totally symmetric. It is straightforward to check 
that the RG equations for the QID systems do not satisfy this criterion 

However, under general linear transformations of the couplings hi[l) = Lijgj^l), 
the coefflcients M^^ transform into a new set of coefficients N^^, which may become 
symmetric. For convenience, we introduce a set of matrices [M{k)]ij = M^^ to 
represent the coefficients. The symmetric criterion for N^^ = N'^^ requires the 
existence of a constant matrix L which satisfies the following constraints (for all fc!), 

M{kY' = [L^L)M{k){L^L)-^, (5-3) 

where superscript T means transpose. In general, there is no guarantee why the 
strongly over-determined constraints would allow a solution for L. In fact, it is a 
nontrivial task to just prove/disprove whether the desired linear transformation L 
exists. Surprisingly, for the QID systems, the hunt for the solution greatly simplifies 
if we formulate the problem in terms of Major ana fermions. The desired transfor- 
mation becomes diagonal, Lij = ri6ij, where rj is a set of rescaling factors and leads 
to the totally symmetric coefficients. 

So the search for the potential is now nailed down to find a set of rescaling factors rj 
in the Majorana representation. In later section, we demonstrate how to construct 
the RG potential explicitly in the Majorana representation. In fact, one can also 
construct the potential V{hi) directly from the RG equations for dopeci^ and half- 
filled^'^' QID systems. Both approaches lead the the same result and the detail work 
will be described elsewhere. 
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It is important to discuss a special set of analytic solution of Eq. (|5-1|) . which is 
closely related to the scaling Ansatz in Eq. (|1-1() . Suppose the initial values of the 
couplings are gi{0) = Gig{0), where ^(0) = U <^ 1 and Gi are order-one constants 
satisfying the non- linear algebraic constraint, 

Gi = M^GjGk. (5-5) 

It is straightforward to show that the ratios between couplings remain the same and 
the complicated equations reduce to single one, 

For repulsive interaction U > 0, the above equation can be solved easily g{l) = 
~ 0) where the divergent length scale Id = 1/U. Note that this implies the 
ratios of different couplings remain fixed in the RG flows, 

9i{l) = (5-7) 

These special analytic solutions are referred as "fixed rays" because the ratios of 
the renormalized couplings remain fixed along the flows. One immediately notices 
that these special set of solutions are nothing but the peculiar Ansatz found in the 
numerics. As explained in the introduction, if the RG potential exists, these fixed 
rays are the asymptotes of the "valleys/ridges" of the potential profile and capture 
the ultimate fate of RG flows completely. 

After familiarizing the readers with the one-loop RG equations, we now proceed 
to construct the RG potential explicitly for QID systems. It turns out that it is 
convenient to represent all interacting vertices in Eqs. H2-27() and (|2-28() by Majorana 
fermions. 

5.1. Majorana Representation 

Now we switch to the Majorana fermion basis and construct the RG potential 
explicitly. Without interactions, the band structure in low-energy limits is described 
by Nf flavors of Dirac fermions with different velocities in general. Each flavor of 
Dirac fermions can be decomposed into two Majorana fermions. Combined with spin 
degeneracy, the 4:Nf flavors of Majorana fermions are described by the Hamiltonian 
density 

Wo = ilRai-'iVadx)'nRa + 'nLa{iVadx)r]La, (S'S) 

where Va denotes the Fermi velocity for each flavor. 

In general, a single vertex involves four different Fermi points, which generally 
would have four different velocities. However, while the vertices in Eq. ()2-27() and 
(|2-28jl may involve four different bands, there are at most two different velocities 
associated with each vertex because Vi = v--, i.e. the velocities in each vertex always 
appear pairwise. This seemingly useless feature turns out to be strong enough to 
guarantee the existence of the RG potential when the Hamitonian is re-expressed in 
terms of Majorana fermions. The momentum conservation in weak coupling somehow 
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V V 

b b' a' b' 
F -vertex B- vertex 

Fig. 1. Forward and backward vertices. Bold lines represent right-moving Majorana fermions while 
dashed lines stand for the left-moving ones. 



gives rise to the interesting constraint that the Fermi velocities must equal pairwise! 
Therefore, the interacting Hamiltonian in terms of the Majorana fermions take the 
form, 

Hint = F{a, a'; 6, h')r]^^r]^^,r]^^r]^y 

+ B{a,b;a,b')r]j,^r]^^7]^^,ri^^,, (5-9) 

where summations over allowed indices are implied. We emphasize that the allowed 
interactions might involve four different Fermi points labeled by a,a',b,b', but only 
two different velocities Va = Va' and Vh = appear in a single vertex. By direct 
comparison, it is clear that the F— vertex include fij and Wij and the i?— vertex cover 
Cij , Sij and Uij . 

There are also other kinds of interactions allowed by momentum conservation, 

He = Ria, a'; b, b')r]j,^r]^^,ri^,r]^^, 

+ Ha, a', b, &')^La^L„'??L6^L5'- (S'lO) 

The scaling dimensions of these vertices are (Z\jj,Z\£) = (2,0), (0,2), while the ver- 
tices in Eq. (|5-9|) have scaling dimensions (1,1). Since the renormalization comes 
from loop integrations, only vertices with scaling dimensions differed by (n, n), where 
n is an integer, would renormalize each other. As a result, the chiral vertices in 
Eq. H5-10p remain marginal and only renormalize the corresponding Fermi veloci- 
ties. Since the corrections only show up at two-loop order, we would ignore their 
contribution here. The pairwise-equal Fermi velocities in Eq. (|5-9() make the classi- 
fication of all vertices fairly simple as shown in Fig. ^ The names come from the 
fact that the forward-type {F-) vertices include the usual forward scatterings while 
the backward-type {B-) vertices include the backward scatterings. 

To obtain the flow equations, we need to integrate out fluctuations at shorter 
length scale successively. The most convenient approach is by the operator product 
expansions we introduced in previous sections. To one-loop order, the renormaliza- 
tion of the bare couplings come from four types of diagrams FF — > F, FB — > B, 
BB — > F and BB B, shown in Fig. [3 Let us start with the first type of loop 
diagrams, FF F in Fig. I2ta). The OPE of Majorana fermions can be computed 
straightforwardly and the mode elimination leads to the renormalized Hamiltonian 
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Fig. 2. Four different diagrams to the one-loop order. (A) FF ^ F (B) FB ^ B (C)BB F (D) 
BB — » B. Notice that, only in the fourth diagram, there are three velocities involved while only 
two velocities are involved in all other diagrams. 



density, 



= F(a, a"; b, b")F(a", a'; b" , b') ''^'^y^'''^>'\^^' dl 

T^iya + Vb) 

= dF{a, a'; b, b')rj^^r]^jr]^,r]^y , 



(5-11) 



where dl = ln6 is the logarithmic length scale. Here we have used the fact that 
Va = Va" = Va' and Vb = Vb" = Vb'- The factor l/TT{va + Vb) arises from the product of 
propagators with opposite chiralities and different velocities. Introducing a simple 
rescaling of the original couplings according to their associated velocities 

F{a, a; b, b') = ^^=F{a, a'; b, b'), (5-12) 



the RG equation is 

dF{a,a';b,b') 
dl 



27:y/VaVb 



^abF{a,a";b, b")F{a" ,a';b" ,b') + ... 



(5-13) 



where ^ab = "^^/VaVb/iva + Vb). Repeating similar calculations, the RG equation for 
F{a, a"; b, b") contains a term 70^^(0, a'; b, a')F{a" , a'; b" , b) and similar result for the 
renormalization of F{a" ,a'; b" , b'). Therefore, the flow equations can be derived from 
a potential. 



ViF) = -jabFia,a';b,b')Fia,a";b,b")F{a",a';b"b'). 



(5-14) 



The RG potential for the {FBB) and (BBB) cases shown in Fig.|2i;b-c) and^d) 
can be constructed in a similar fashion. Finally, combining all contributions together, 
the weak-coupling RG flows for QID systems are described by the potential 



V{F, B) = -B{a, b; a, b')B{b, c; b' , c)B{c, a; c , a 

— ^abB{a, b'; a, b)B{a" , b'; a , b")F{a, a"; b, b") 

- jabFia, a'; b, b')F{a, a"; b, b")F{a",a'; b"b'), 



(5-15) 
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Summations over all allowed indices are again implied. The merits to use Majorana 
representation enables us to construct the RG potential explicitly. Compared with 
previous examples for the half-filled two-leg ladder and the doped A^-leg system at 
generic filling, one can easily see that the RG potentials in Eqs. (|4-14|) and (|4-17() 
agree with the general form in Eq. (|5-15|) . Note that if one starts from a 'wrong' basis, 
it is far from trivial to realize the fact that the non-linear flows can be derived from a 
single potential. However, in the Majorana representation, the linear transformation 
to the potential basis is diagonal Lij = ViSij. After appropriate rescaling of couplings, 
the explicit form of the RG potential is derived. We have checked for the doped and 
half-filled QID systems and found all potentials agree with Eq. (|5-15|) . 

There is one loose end about the rescaling factors. For most physical systems, 
the rescaling factor is slightly more complicated than {2'iry/VaVh)~^ ■ the subtlety 
arises from the degeneracies of the couplings imposed by physical symmetries. This 
is best illustrated by the following simple example. Consider the RG equations for 
three couplings gi, where i = 1,2,3, 

Since |eijfc| is totally symmetric, the corresponding RG potential is V(^g) = 515253- 
Suppose the system has some symmetry, such as U{1) symmetry for charge conser- 
vation, and the couplings are degenerate 52 =53- The RG equations are simplified, 

dgi 2 dg2 

^ = 52, ^=9192. (5-17) 

It is straightforward to show that we need to perform a rescaling transformation, 
{hi, /i2) = (51, \/252) to obtain the potential V{h) = /ii(/i2)^/2. In fact, for couplings 
with re-fold degeneracy, an additional rescaling factor ^/n is necessary to bring them 
into the potential basis. Therefore, the total rescaling factor is 



n = ^\f^. (5-18) 

where rej is the degeneracy number of the coupling gi, with Fermi velocities Va and 

Vb- 

So far, we have shown the existence of the RG potential by explicit construction, 
which proves the the widely used Ansatz in Eq. (jl-lf) . In fact, the asymptotes of 
the RG flows are governed by the special set of fixed-ray solutions. This in turns 
explains the simplicity of the phase diagram, even though the RG equations are 
rather complicated. 

The absence of exotic fates of the RG flows is also found in earlier work on 2D 
melting theorjISSi qj- Kondo related problems However, the simplicity of RG 
flows in these systems is not quite the same as described here. Since the number 
of marginal couplings in these problems is few, analytic solution of the flows often 
show that it is possible to define some conserved quantity associated with each flow 
lines. This is where the simplicity comes from. On the other hand, we have also 
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looked into these well-known flows to check whether they can be derived from a 
single potential. It is not too surprising that this is indeed the case because the 
requirement of potential flows loose up quite a bit when the number of couplings is 
small. 

With the help of non-perturbative Abelian bosonization, it is not essentially im- 
portant whether the RG flows can be cast into potential form beyond one-loop order. 
However, it remains an interesting and open question at this moment. Note that 
the coefficients of the one-loop RG equations are unique, protected by the leading 
logarithmic divergences. The next order calculations bring in lots of complications 
and subtleties, including the non- universal coefficients in the RG equations, velocity 
renormalization and so on. It is not clear at this moment whether it is even sensible 
to pursue the RG potential beyond one-loop order. 

Another interesting issue concerns the connection between the potential V{hi) 
and the Zamolodchkov's c-function C{gi) of (l-l-l)-dimensional systems with Lorentz 
and translational symmetries.!^ A generic QID system we studied here has neither 
Lorentz invariance (due to different Fermi velocities) nor translation symmetry (due 
to Umklapp processes). While both V and C are non-increasing along the RG 
flows, the exact relation between them remains unclear at this point. We emphasize 
that the existence of a non-decreasing function C along RG flows only implies that 
dC/dl = (dC/dgi) ■ {dgi/dl) < and is not strong enough to show that the flows 
can be derived from a potential. Thus, the potential flows are closely related to the 
c-theorem but they are not equivalent in general. In addition, we do not know any 
easy generalization of c-theorem that does not rely on the Lorentz and translational 
symmetries. However, one can easily check that in the special limiting case where 
Lorentz and translation symmetries are restored, the C-function indeed coincides 
with the potential we find. This indicates that there may be a general form of 
c-theorem waiting to be discovered. 

In summary of this section, we have shown that the RG transformation for QID 
systems in weak coupling is described by potential flows. Therefore, neither chaotic 
behaviors nor exotic limit cycles could occur. The different Fermi velocities and 
the degeneracies imposed by physical symmetries give rise to non-trivial rescaling 
factors, which hinder this beautiful structure behind the RG transformation. The 
explicit form of the potential is obtained after appropriate rescaling of the couplings 
in Majorana basis. 

§6. Anomalous Scaling 

Now that we have established the existence of the RG potential, we would like 
to discuss its physical consequences. Since the notion of RG potential is still novel to 
the community of low-dimensional correlated physics, we concentrate on the simplest 
example, the half-filled two-leg ladder, in this section. However, one should keep in 
mind that most of the results discussed here can be applied to more general ladder 
systems. 

For the two-leg ladder at half filling, the number of possible interactions is greatly 
reduced to nine in weak coupling. Within the one-loop RG calculations, these cou- 
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plings gi are described by a set of coupled first-order differential equations as in 
Eq. (|5-1|) . For the two- leg ladder at half filling, Mj^ are 9x9 constant matrices, 
as implicitly contained in the RG equations in Appendix B. For a generic interact- 
ing Hamiltonian, the bare values for these nine couplings can be straightforwardly 
determined. However, it is generally very difficult to obtain the solution for these 
coupled flow equations in analytical form. 

Simple analytical solutions emerge if the interactions are chosen in a specific 
way. These special solutions are later referred as "symmetric rays". Suppose the 
bare couplings of the specific interacting Hamiltonian are gi{0) = Gig{0), where 
5(0) = (f//t) <C 1 is small while Gi are order one constants which satisfy the algebraic 
constraint, 

= MjfcG.Gfe. (6-1) 

It is straightforward to show that the ratios between couplings remain the same 
and the nine complicated equations reduce to single one, g = g'^ I The solution is 
g{l) = l/{ld — I), where the divergent length scale Ij, = {t/U). Of course, one should 
keep in mind that the solution g{l) is only valid when it does not flow out of the 
weak coupling regime. These special Hamiltonians, whose couplings are described 
by these symmetric rays, turn out to be S0(8) symmetric. For the two-leg ladder, 
four different phases,»^''^' named as D-Mott, S-Mott, CDW and SF, are of the central 
concerns. 

It was shown previously that the two-leg ladder in weak coupling always scales 
into one of the four different symmetric phases.^ However, this numerical approach 
was criticized that the SO (8) symmetric phases might have instabilities which happen 
not to be tackled by the limited types of interactions considered in the numerical 
study. To make up the fissure, a complete stability check near the S0(8) symmetric 
rays is desirable. 

6.1. Stability Analysis 

To describe the RG flows in the vicinity of the S0(8) symmetric rays, it is 
sufficient to consider the linearized version of Eq. ()5-l() . For a generic interaction, 
the couplings are separated into symmetric and asymmetric parts, gi{l) = Gig{l) + 
Agi{l). In the vicinity of the symmetric rays, the deviations are small, Agi{l) <^ g{l)- 
Keeping the leading order term, the linearized RG equations are 



diAgi) _ Bij 



dl {Id - I) 



^g^, (6-2) 



where Bij = 2A^-fj-k- The matrix Bij can be brought into diagonal form by a linear 
transformation. As a consequence, the RG equations decouple into nine independent 
ones, 

d{5gi) > 



dl {Id - I) 



Sgi, (6-3) 



where 6gi are couplings after the linear transformation and Aj are the eigenvalues 
of the matrix Bij. Although the matrix Bij are different for each S0(8) symmetric 
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rays, the eigenvalues are identically the same 



X^ = 2, 



2 2 2 



3 3 3 



1 

3' 



1 

3' 



1 

'3' 



(6-4) 



This coincidence implies that the results of the stability check only rely on the 
symmetry group but not on the details of the phases 

So far, we have a single equation, g = g^, describing the renormalization along 
the symmetric rays and nine for deviations from the rays as in Eq. 1)6 -3^ . Apparently, 
there must be one redundant equation among them because the original number of 
equations is only nine. It doesn't take long to find out that it corresponds to the 
flow equation with largest eigenvalue A = 2. This corresponds to a trivial case that 
all couplings are shifted along the symmetric ray, i.e. dgi = Gi5g. Since the flow 
along the ray is described hy g = g^, linearization leads to 6g = 2g{l)5g. Since we 
only need eight equations to describe the deviations from the symmetric rays, the 
A = 2 equation is only an artifact and should be ignored. 



(a)X<0 



(b) 0<A<1 



(c) ?L>1 



Fig. 3. The topology of RG flows near the symmetric ray with (a) A < 0, (b) < A < 1 and (c) 
A > 1. It is clear that the coupling is irrelevant for A < and relevant for A > 1. The RG flow 
is more subtle for < A < 1. In this case, although the deviation from the symmetric ray is 
growing, the slope remains the same as the ray. 

The other eight eigenvalues describe how the flow goes once the bare couplings 
are off the ray. Starting from the bare values Sgi{0) = 5Gig{0), where 6Gi <C Gj. 
The solutions of Eq. (|6-3() are 



Sgiil) 



5Gi 



{Id - ly 



(6-5) 



According to conventional classification,'^ for Aj < the deviations diminish under 
RG transformations as shown in Fig. Ola) and thus are classified as irrelevant cou- 
plings. For Aj > 0, any small deviations get enhanced and the flow is pushed away 
from the symmetric rays as in Fig. ISl^b) andEJc). These are classified as relevant 
couplings. The conventional wisdom tells us that there are three relevant couplings 
(Aj = 2/3) and five irrelevant ones (Aj = —1/3). Since a generic interaction in prin- 
ciple could generate asymmetric deviations in all couplings, one might rush to the 
incorrect conclusion that the S0(8) symmetry is not stable. However, the first-glance 
guess is wrong because the conventional classification is based on the perturbative 
analysis near a fixed point while we are dealing with running symmetric rays. A new 
set of rules to identify relevant perturbations is in order. 
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The crucial criterion is whether the deviations Sgi{l) grow larger than the sym- 
metric coupling g{l) = (l^ — Just growing larger than the bare values under 
RG transformations is not qualified as a relevant perturbation. This subtle but 
important difi'erence is best illustrated by calculating gap functions using scaling 
arguments. If the degeneracy of excitation gaps is maintained, the S0(8) symmetry 
is robust and vice versa. 

6.2. Scaling Argument and Anomalous Exponents 

Now we elaborate on the anomalous scaling arisen from the peculiar RG fiows 
near the symmetric rays. Let us remind you that there are several useful properties 
of the RG flows in Eq. (|5-1|) . Suppose all bare couplings are rescaled by a factor of 
U, gi{0) = Ugi{0), it is easy to check that the RG fiows are recovered by a rescaled 
(logarithmic) length scale I = l/U. That is to say, 

~gi{l) = UgiiUl). (6-6) 

It means that, as long as the relative ratios of the bare couplings are fixed, the RG 
fiows for different interaction strength are related in a simple way. 

To gain some physical intuition of the scaling argument, it is illuminating to 
study the simplest one-loop equation 

-dl=-'^ ^'-^^ 
with solution g{l) = 1/(1 — Id), where the divergent length scale Id = —l/g{0). 
The relevance of the marginal coupling g depends on the sign of the bare coupling 
value. For g(0) > 0, Id is negative so that the renormalized coupling g(l) scale 
toward zero as power law. For g{0) < 0, g{l) grows under RG transformation and 
eventually becomes divergent at / = Id- This simple RG equation actually describes 
the well-known BCS instability. For repulsive interaction, corresponding to g{0) > 0, 
the pairing interaction is marginally irrelevant and the Fermi liquid is stable with 
inclusion of electronic correlations. On the other hand, for attractive interaction 
where g{0) < 0, the pairing instability is enhanced and leads to superconducting 
ground state. 

The analytical solution in the previous example also help us to understand the 
peculiar scaling of the gap function generated by marginally relevant couplings. We 
must emphasize that the appearance of a relevant coupling does not necessarily 
mean the formation of some excitation gap in the energy spectrum. However, let us 
concentrate on the case where the relevant coupling does drive the formation of some 
gap. From conventional scaling argument, Z\[g(Z)] = e'z\[(7(0)], we cut off the fiow 
when the relevant coupling reaches one g{lc) = 1 because the one-loop RG equation is 
no longer reliable. It is easy to solve for the cutoff length scale, lc = ld — ^^ l/lffCO)!- 
At the cutoff length scale, since the coupling is order of unity, the gap function is 
also order one. Thus, the scaling argument tells us that 

A[giO)] = e-'^A[g{l,)] ^ Aq x e-VI^^WI, (g.s) 

where Aq is some is some energy scale of the system (like the band width t). Although 
the non-analytical gap function is derived in this extremely simple case, it is expected 
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that the gap function would behave in a similar way for more general quasi-lD 
systems. 

Since the RG equations are only valid in weak coupling, we cut off the RG 
procedure when the symmetric coupling g{lc) = 1- At this cutoff length scale Ic, the 
deviations in Eq. H6-5() are 



Sdiilc) ~ [jj ■ (6-9) 

For Aj > 1, the deviations are larger than the symmetric coupling and should be 
classified as "relevant" . As long as Aj < 1, the deviations at the cutoff length scale are 
still vanishingly small and should be viewed as "irrelevant". This new classification 
is different from the conventional one for < Aj < 1 . We would see clearly soon that 
the new classification is appropriate for stability check near a running symmetric 
ray. The grey regime < A, < 1 between the new and conventional rules gives rise 
to anomalous scaling exponents. 

Under RG transformations, the gap functions scale like, 

A [g{0),6gim = e~'=A [l,Sgi{lc)] , (6-10) 

where 6gi{lc) are at most of order (U/ty^^. The key point is that, although some 
perturbations are enhanced to order (U/t)^^^, which is larger than the bare order U jt 
values, they are still small at the cutoff length scale. The effective Hamiltonian can 
be separated into two parts H = Hq + 6H - the SO (8) symmetric and asymmetric 
parts. Since the asymmetric part is of order {U /t)^/^ , the changes of the gaps would 
be of the same order by standard perturbation theory. Without the deviations, the 
S0(8) symmetry guarantees the exact degeneracy of all gaps, i.e. /1[1,0] = A. The 
presence of perturbations modifies the gap functions, 



Ai[l,5gi{Q]= A 



1 + Ci ( y V + 



(6-11) 



It is clear that, in the weak coupling limit U/t 0, the degeneracy of all gaps 
is recovered. It implies that the S0(8) symmetry is indeed robust under generic 
perturbations. 

The anomalous scaling exponent 1/3 is clearly seen in the ratios between charge, 
spin and single particle gaps (see Fig. [l]) , 

^«>l+c«(^)'. (6.12) 

This simple exponent is rather fascinating because the RG equations we started from 
are very messy. But, despite of which S0(8) symmetric phases the system flows into, 
the exponent of the gap function corrections is universally equal to 1 — Aj = 1/3! 

It is quite exciting to learn that the existence of RG potential, not only explain 
the simplicity of the phase diagrams from the complicated RG equations, it also 
predicts the emergence of the enhanced symmetry in QID ladders. Furthermore, 
the instability check in the vicinity of these symmetric rays gives anomalous scaling 
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Fig. 4. Ratios of charge, spin and single particle gaps plotted versus the interaction strength U /t. 
The ratios approach unity in the asymptotic U/t ^ limit with anomalous corrections of order 



with ratios numbers associated with the emergent symmetry. While we only worked 
out the details for the half-filled two-leg ladder in this section, it is expected that 
the phenomena should be rather general. There is still plenty of room to investigate 
the fun physics associated with the RG potential in QID ladders. 

§7. Discussions and Conclusions 

In this article, we reviewed how to write down the effective action for QID 
systems in weak coupling. In contrast to conventional 5f-ology, we classify and express 
all possible four-fermion interactions by SU(2) bilinear currents. With the help of 
OPE's, the derivation of RG equations becomes equivalent to obtain the coefficients 
of algebra among the currents. To our best knowledge, this is the first time that the 
whole set of RG equations for generic QID systems is written down explicitly. 

While numerical elaboration of the RG equations is important to reveal physical 
properties in the ladder materials, we focus on the hidden structure behind the 
scene and hope to answer two fundamental questions: why the phase diagram is 
so simple despite of the complexity of the non-linear RG equations and why the 
enhanced symmetry emerges after integrating out the short-range fluctuations. Both 
puzzles can be answered by the non-trivial existence of the so-called RG potential 
which allows a simple geometric interpretation of the flows. The ultimate destinies, 
after successively integrating out fluctuations at shorter length scale, are nothing 
but the well- isolated asymptotic "valleys/ridges" with specific ratios between the 
renormalized couplings. The isolation of the asymptotes explains the simple and 
robust structure of the ground-state phase diagram and the specific ratios give rise 
to the enhanced symmetry emerging in low-energy limit. We first demonstrated how 
the RG potential can be constructed explicitly in the half-filled two-leg ladder and 
also the doped A^-leg ladder. Then, we provide a general proof for its existence for 
generic QID systems and also pin down the functional form of the potential. 

There are many open questions left behind the work here. For instance, it 
is interesting to filling in the numerical details for specific ladder materials. But, 
maybe one of the most interesting questions is its applicability to strongly correlated 
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systems in two dimensions. First of all, one may naively expect QID systems with 
large N should share a similar ground state for a truly two-dimensional system. 
Unfortunately, if the Fermi surface is not nested, the instability in the RG equations 
fade away and the ground state is the normal Fermi liquid. The failure to deliver 
an exotic ground state in the presence of featureless Fermi surface may indicate two 
possibile scenarios. 

The first scenario is that all interesting ground states in two dimensions only 
occur in strong coupling and the RG analysis in weak coupling provides null informa- 
tion about them. If so, one should abandon the weak-coupling approach completely 
and search for alternative routes. On the other hand, the second scenario to explain 
the failure is tied up with the shape of Fermi surface. It is possible that the exotic 
phases in two dimensions are driven by the enhanced couplings due to nested Fermi 
surface. Therefore, by cutting up the Fermi surface appropriately, the RG technique 
developed here can be applied to determine the ground states in two dimensions. 
Even if the instabilities towards some exotic phases can be spotted in weak-coupling 
RG, it doesn't mean that we understand or know how to appropriately describe 
the physical properties which axe essentially described by the fixed points in strong 
coupling. However, the powerful technique of bosonization may help us out here to 
sketch the ground state qualitatively. It is not clear at this stage which scenario is 
correct and more numerical work needs to be done to clear up the relevance between 
the QID and the true 2D physics. If the truth favors the second scenario, the RG 
potential would give us more confidence for what is found in weak-coupling analysis. 
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Following similar calculations presented in the text, one can obtain the full set 
of the OPE's. For interested readers who would like to derive the identities on their 
own, you will find the following identities useful. 
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Appendix A 

Operator Product Expansions 



(A-l) 



In below, we list the complete rules for the OPE of right-moving currents. 




(A-2) 



32 Ming-Shyang Chang, Wei Chen and Hsiu-Hau Lin 

J5(x,r)J,K0,0) ~ ||J^ - ^J,«, (A.3) 



+^^'^'^1^^^ + ^^^,^ (A-4) 



lf^{x,r)Um - ^^J!i + ^ - ^^Ja - ^Jfk, (A-6) 

l2;;'^^'"2;;'^^'="2;"'^^' + 2i"'^^^/' 



hj{x,r)Jki{0,0) ~ ^Iji + ^lu, (A-8) 

h,{x,T)Jum - ^7r.^!i + ^7^11 (A-9) 

/^.(x,r)JH(0,0) ~ -^/« + (A-IO) 
There are other OPE's, for example, ijjJki, can be obtained by taking hermitian con- 



jugate in Eq. (|A-8|) on both sides and sending Zj — > —Zi at the same time. Replacing 
Zi — > z* in the above equations gives the OPE's for left moving currents. 

Appendix B 

RG Equations 



Since the four-fermion interactions in QID system are all marginal, the renor- 
malized couplings are described by coupled non-linear equations. In addition to the 
usual forward and Cooper scatterings, denoted as f[j^ and c^j"^, the nested Fermi 
surface also host the extra momentum-conserving vertices s^j*^ and the umklapp 
interactions Uif, w^'" . After integrating out fluctuations at short-length scale, the 
complete RG equations for the QID system at half filling are shown in the followings: 
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k ^ ^ 

= ~ ^^^\j^ik^kj + CikCkj + "^Clk^kj^ + -^ic^jKj + Ci^i'^fj - '^Cijhfj) 

/^=i[(4)'+3(4)^; 

+ E^.i^{(4)^ + 3«.)^ + Hkf + K<kf] 

+ ^fi<? - ^[(^g)' + 3.^.)'] + ^[(^g)' + 3(«|)^] + + 3K)^], 

(B-3) 

Jij \Jtj) ^ 2 *J '-^ 2 

+ 2[^^^^ - - + + .^.^ + [s^] (B.4) 

k ^ ^ 

+ \K< + + 34i^- + 4< + 3^5c^), (B-5) 

^ij = E ~4^|'^ifc^fci + ^ik^kj ~ '^^Ik^kj + "ifc^fej + '^ik'^kj ~ '^'^ik'^kj^ 

+ \^^^ + 4< - 24^^ + ^^0 + ^4^" + - 24/^^)' 

(B-6) 
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+ 



u 



(B-7) 



— X/ ^^"i ^''^ik^j '^^fc^fcj + '"'ik^kj + '^ifc^fcj ~ '^4k''^kj ~ '^''Ak^j) 



k 



(B-8) 



1 



(B-9) 



Notice all couplings in the RG equations are rescaled with a velocity factor to simplify 
the coefficients, gij = gij/'iT{vi + Vj). The dimensionless constant aij^h inside the 
summations depends on the Fermi velocities, 

Note that oiij^k = 1 when Vi = or vj = v^- Finally, several short-hand notations 
are defined, 

hij = 2fij + SijCii, (B-12) 
u^j = Uii + Ujj, (B-13) 
4^^ii + %- + 2pAi' (B-14) 
where A = {c, /, s} and p = 0, ±1. 
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